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Abstract 

The effects of bulk viscosity on the elliptic flow V2 are studied using realistic equation of state 
and realistic transport coefficients. We find that the bulk viscosity acts in a non trivial manner on 
V2- At low pr , the reduction of x>2 is even more effective compared to the case of shear viscosity, 
whereas at high pj, an enhancement of v>2 compared to the ideal case is observed. We argue that 
this is caused by the competition of the critical behavior of the equation of state and the transport 
coefficients. 



1. Introduction 

The effect of shear viscosity in heavy ion collisions has been studied by several authors Hj] 
and was found to affect considerably the elliptic flow observable \>z. On the other hand, the bulk 
viscosity has always been neglected and its effects are mostly unknown. The shear viscosity 
acts as the resistance against the deformation of a fluid element. In this sense, it is natural to 
expect that va, which characterizes the spatial anisotropy of the dynamics, is affected by the shear 
viscosity. On the other hand, the bulk viscosity acts only against the expansion or compression 
of the fluid. However, this does not mean that the bulk viscosity has negligible effects on Vz, 
particularly because the expansion and compression of the produced matter in the collisions is 
also anisotropic. 

Therefore, the behavior of the bulk viscosity coefficient, near the QCD phase transition 
should also be carefully studied. Recently, lattice QCD simulations J2l suggest that f is strongly 
enhanced near the critical temperature T c , A similar enhancement will appear even below T c , as 
is shown in yfl, where £ in the hadronic phase increases also quickly towards T c , Considering 
the fact that the fluid dynamics slows down near T c due to the small sound velocity c s , the critical 
behavior of f should have a crucial role on the expansion of the matter. 

In this work , we investigate the effect of bulk viscosity on V2 by incorporating the critical 
behavior of the equation of state (EoS) and £ within a 2+1 dimensional hydrodynamic modeling. 
On the other hand, the dissipative corrections to the one-particle distribution function at freeze 
out is not considered, because Grad's method is not applicable for the freeze out values of bulk 
viscosity obtained in our calculations |4]. 
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2. Relativistic Dissipative Hydrodynamics 



The formulation of relativistic dissipative hydrodynamics is not trivial. The naive relativistic 
generalization of the Navier-Stokes equation violates causality and leads to intrinsic instabilities 
making the theory not applicable. In this work, we use the memory function method |3,|g], which 
is formulated in such a way that the magnitude of the bulk viscosity has a lower bound, and it is 
possible to implement stable numerical simulations even for ultra-relativistic initial conditions. 

We consider the case of vanishing baryon chemical potential where only the conservation of 
the energy and momentum is required. For a general metric g^ v , the continuity equation for the 
energy-momentum tensor is given by 

^^(V^J + ry^o, (i) 

where V. is the Christoffel symbol and g is the determinant of g^y. We use the Landau definition 
for the local rest frame and assume, as usual, that the thermodynamic relations are valid in this 
frame. We further ignore the shear viscosity. Then the energy-momentum tensor is expressed as 
T 1 ™ — (s + p + n) m /j m v - (p + IT) g^ v , where s, p, u M and II are the energy density, pressure, four 
velocity and bulk viscosity, respectively. The equation for the entropy production is, 

-La, ( v^) = - 5 -La, ( v=^) - (2) 

where s is the entropy density. From this, we can define the thermodynamic force F associated 
with the bulk viscosity as ■yj—g~ l d tl y yf^gurj. It is also useful to introduce the specific volume <x 

defined by the conservation law yf-g^d^ ( -^gcmr \ — 0. 

In the memory function method, an irreversible current / is induced by the corresponding 
thermodynamic force as 



m 

ct(t) 
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Thus, the equation of the bulk viscosity is obtained by setting J — II and F = yf^g 1 ( yf^gurj. 
This is re-expressed in a differential form as follows, 

TKU"d„ (-) + - = ( \Fg« f ') , (4) 

Wl a- a- V=g v ' 

where tr is the relaxation time and £ is the bulk viscosity coefficient. 

In the following calculations, we consider the hyperbolic coordinate system = (t, r^, 77), 
where t = \/t 2 - z 2 , tj - \ tanh(^jf), r T = (x,y) and yf^g = t. 

To solve numerically the relativistic hydrodynamic equations ([T]), (f2| and ©, we use the 
Smoothed Particle Hydrodynamic method The initial condition used is a parametriza- 

tion of the usual Glauber model. The form of this initial condition and the choice of parameters 
can be found in [4]. The impact parameter b is taken to be 7 fm. For the EoS, we consider a 
smooth interpolation between the lattice QCD results Ht] and the hadron resonance gas, which is 
shown in Fig. Q] 

For the bulk viscosity coefficient, we use the result from lattice QCD J^l for the QGP phase 
and the result from |3[ for the hadron phase. Then, the ratio g/s exhibits a maximum at T c and 
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starts to decrease almost exponentially as the system departs from T c . We express these features 
of f/i as 

( A l x 2 +A 2 x-A i (0.99570 > 70 > 1.05T C ) 

- = < A] exp(-(jc- l)/<ri) + ^ 2 exp(-(x - l)/cr 2 ) + 0.001 (J > 1.05r c ) (5) 
s { A 3 exp((jc - l)/cr 3 ) + A 4 exp((jc - l)/cr 4 ) + 0.03 (T < 0.99570 c ), 

where x = T/Tc- The fitted parameters are A\ = A-$ = 0.9, /I2 = 0.25, /I4 = 0.22, en = IO0-3 
0.02S. 0-4 = 0.022, A\ = -13.77, A 2 = 27.55 and A 3 = 13.45. For the relaxation time, we use 
the following parameterization t r - bn(/(e + p) @]. In this work, we use b n = 6 J3l. 



3. Numerical results 

We calculate the elliptic flow parameter, V2 = (cos(20)), for pions as a function of the trans- 
verse momentum, using the Cooper-Frye method with a freeze out temperature 130 MeV fl. 
In Fig. |2] v'2 is shown for viscous fluids (dashed, dash-dotted and dotted lines) and ideal fluid 
(solid line). Dashed and dash-dotted lines correspond to the constant g/s cases and are shown 
just for the sake of comparison, whereas the dotted line represents the £/s with the QCD phase 
transition. 

For small transverse momenta Pj, V2 is suppressed by the effect of bulk viscosity, similar 
to the case of shear viscosity. However, at high Pr, V2 starts to recover and even surpasses the 
values obtained from ideal hydrodynamics. This effect is enhanced when £/s exhibits a critical 
behavior. This behavior at high Pj is not observed for the case of shear viscosity and we further 
confirmed that this is not observed when an ideal gas of massless quarks EoS is used 041]. Thus, 
this effect should be attributed to the presence of the QCD phase transition, both in the EoS and 
in the bulk viscosity. 

This behavior can be understood as follows. For the viscous case, the radial dependence 
of the fluid flow is not monotonic in the region where the velocity of sound is minimum. This 
is because the flow of the internal matter, which has a higher temperature, tends to catch up 

3 




p T (GeV) 

Figure 2: v>2 as a function of pj. The figure displays results for different choices of f/i, The solid line corresponds to 
the ideal result, the dashed and dash-dotted lines correspond respectively to (/s = 0.08 and 0.2, and the dotted line is the 
result for f/j with critical behavior. 

the foregoing fluid elements, generating an inhomogeneous velocity field in the radial direction. 
Then, the bulk viscosity n starts to increase in this region, heating the matter. If this happens, 
the gradient of PI becomes dominant in the acceleration of the fluid compared to the pressure 
gradient (note that the acceleration is given by -V (p + IT) ) since the gradient of the pressure 
is proportional to the sound velocity. Such a mechanism works more efficiently in the direction 
where the initial acceleration is large, and in consequence, in the direction to increase the elliptic 
flow. Furthermore, this effect becomes effective only when significant collective flow is formed 
near T c . Thus the recovery of V2 by this mechanism is expected for higher momentum particles. 

We acknowledge illuminating discussions with T. Hirano and A. Monnai. This work has 
been supported by CNPq, FAPERJ, CAPES, PRONEX and the Helmholtz International Cen- 
ter for FAIR within the framework of the LOEWE program (Landesoffensive zur Entwicklung 
Wissenschaftlich- Okonomischer Exzellenz) launched by the State of Hesse. 

References 

[1] See, for example, the contribution from P. Romatschke in this volume. 

[2] F. Karsch, D. Kharzeev and K. Tuchin, Phys. Lett. B663, 217-221 (2008). 

[3] J. Noronha-Hostler, J. Noronha and C. Greiner, arXiv:0 8l 1.1 57 11(20 08) 

[4] G. S. Denicol, T. Kodama, T. Koide and Ph. Mota, arXiv:0903.3595 (2009). 

[5] T. Koide, G. Denicol, Ph. Mota and T. Kodama, Phys. Rev. C75, 034909 (2007). 

[6] G. S. Denicol, T. Kodama, T. Koide and Ph. Mota, J. Phys. G35, 115102 (2008). 

[7] G. S. Denicol, T. Kodama, T. Koide and Ph. Mota, Phys. Rev. C78, 034901 (2008). 

[8] Y. Aoki et al (WHOT-QCD Collaboration), arXiv:08 10.4742 (2008). 



4 



